Phase Transition in a Self— Gravitating Planar Gas 
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We consider a gas of Newtonian self-gravitating particles in two-dimensional space, finding a phase 
transition, with a high temperature homogeneous phase and a low temperature clumped one. We 
argue that the system is described in terms of a gas with fractal behaviour. 



I. INTRODUCTION 

The statistical behaviour of systems interacting via 
classical (Newtonian) gravitational forces is rather pe- 
culiar as compared with other statistical systems, such 
as neutral gases and plasmas. The central distinguishing 
feature of the former is the fact that there is no shield- 
ing of the long range gravitational force, while Debye- 
screening rules the long distance behaviour of the electric 
Coulomb force. In lower dimensions such a characteristic 
of the gravitational force is dramatic, due to the rising of 
the two body potential in Newtonian gravity. This leads 
to conceptual problems, related to the extensive nature 
of energy, as required in statistical mechanics. 

In general these statistical problems are not percept- 
able to exact treatment, and only numerical results are 
known. While classical gravitating systems in one space 
dimension are quite pathological, in two space dimen- 
sions they become manageable, although the most realis- 
tic three dimensional case is still beyond reach. In partic- 
ular, in the latter case the phase space volume diverges 
and one has to use a short distance cutoff. 

In two-dimensional space, the thermodynamic func- 
tions are analytically computable. The thermodynamical 
properties of a gravitational gas has been studied in de- 
tail both from a theoretical jy as well as from a numerical 
@-§ point of view. In general, it has been shown 1gIJ7|] 
that a classical gas of gravitating particles, in the grand 
canonical formalism is related to a field theory described 
by a Liouville Lagrangian. In a two-dimensional space 
Liouville theory is well known Q and several correlation 
functions can be exactly computed |p|-pl| . Indeed, in this 
particular case the theory is conformally invariant, and 
correlators can be computed in terms of known functions 
of mathematical physics, once properly regularized Q . 

At low temperatures the gravitational force is strong, 
that is, the potential energy is big as compared to the ki- 
netic energy, giving rise to a collapsing phase, identified 
by the presence of a single cluster of particles floating 
in a diluted homogeneous background. At high energy 
a homogeneous phase is recovered and the cluster disap- 



pears. Using the microcanonical ensemble one can show 
that in the transition region the system is characterized 
by a negative specific heat (the corresponding instability 
is of extreme relevance for astrophysics, p"3|). 

Hertel and Thirring Q have shown that the canonical 
and microcanonical ensembles are not equivalent in the 
proximity of the transition. This thermodynamic incon- 
sistency has been solved in [[I], and these results have 
been successfully confirmed by numerical investigations 
on self gravitating non-singular systems with short range 
interactions 

Recently, using numerical calculations, a long range 
attractive potential has been considered, as constructed 
taking the first few terms of the Fourier expansion of the 
logarithmic potential. As a result of the simulations, the 
system turns out to exhibit a transition from a collapsed 
phase at low temperature, to a homogeneous phase ||. 

Here we consider a two-dimensional gas of classically 
interacting particles via the Newtonian logarithmic po- 
tential, computing its thermodynamical properties. Us- 
ing results known in the framework of two-dimensional 
conformally invariant euclidian field theory, the complex 
integrals may be computed in closed form, and we find a 
sofisticated structure of poles and zeros in terms of the 
temperature. Furthermore, we introduce an order pa- 
rameter, corresponding to the expectation value of the 
square of the two-body distance. We find that at least at 
one of the singularities of the partition function the aver- 
age distance vanishes, signalizing a clumped phase. We 
are thus able to investigate the collapse of the gravitating 
system. The phase transition point is given in terms of 
the mass of the particles (m) and the gravitational con- 
stant G by T c = j^NGm 2 . We find the exponents of the 
phase transition. 

We furthermore discuss the formulation of gravity in 
terms of the grand-canonical partition function and the 
related Liouville theory as proposed by de Vega et al [|| , 
outlining the consequent relations. We find, using the 
same type of complex integrals, a critical temperature in 
this formulation. 
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II. CANONICAL PARTITION FUNCTION 

We consider a gas of nonrelativistic particles with mass 
m interacting through Newtonian gravity at temperature 
T. Here, the number of particles is fixed and we shall 
work with the canonical ensemble. The canonical parti- 
tion function of the system can be written as 



N 



[d 2 n][d 2 Pl 



(1) 



where [d 2 ri] = HfL-^cPri, and the iV-particle Hamiltonian 
Hn is obtained adding the Newtonian potential to the 
usual kinetic term, that is, 



H 



N 



N y 

^ 2m 

i=l 



l^Gm 2 log| 



(2) 



The potential term must be reqularized, since in gen- 
eral there are divergences in either the infrared or in the 
ultraviolet domain. This is done implicitely definning 
the two-dimensional integrals in the complex plane, as 
done in reference [|l2| (see also ||). We thus consider 
the two-dimensional variables as one single complex vari- 
able in the complex plane, denoting the procedure by 
the use of latin letters from the end of the alphabet, i.e., 
d 2 ri — > d 2 Wi. As it turns out, there is an underlining 
conformal invariance in the present problem, easily seen 
using the complex variables, where is is characterized by 
a SL(2,C) invariance of the Hamiltonian. Such an in- 
variance is valid in the infinite volume limit. For finite 
volume one should take into account boundary effects, 
and in that case it would not be possible to treat the 
problem exactly. Fortunately, due to the fact that gravity 
on a plane constitutes a long range attractive force, the 
system is naturally bound to a region defined in terms 
of the typical bound state lenght. We thus expect the 
procedure to be exact for not too large values of the tem- 
perature. Therefore, well above the critical point depar- 
ture might be expected from our results. Nevertheless, 
the critical temperature is expected to be correct on the 
above grounds. 

Using that invariance, we set wn+i = 0, wn+2 = 1 
and wn+3 — oo (whose integration would lead to an infi- 
nite overall factor in the partition function which can be 
discarded in the computation of physical quantities) and 
find 



z = (2TTm) N + 3 J d 2 Wl \w t \ 2a \i-w t \ 2 ^u l<J \ 



4/< 



(3) 



where a = j3 = — ^Sr and 2p = a. We keep the use 
of the a, (3 and 7 parameters in order to facilitate com- 
parison with the corresponding integrals in the literature 
H[l2]]. Using such well known formulae, the partition 
function Z can be written in terms of p as 



Z^T(N + l)[A(l-p)] N Y[A(zp)x 

1 

JV-l 

[] [A(l + (i + 2)p)] 2 A(-l - (N + 3 + i)p) (4) 



where A (a;) 



This function has poles at the 



r(x) 
r(i-x)- 

points x = —n = 0, —1, —2, • • • and zeros for x = n = 
1, 2, • • •. This is an exact result, as shown in e.g. [|l2|,[s||. 
However, in order to obtain the singularity structure 
(poles and zeros) in further related integrals, it is worth- 
while to classify them. They arise from the following kind 
of rather naive argumentation. 

If one of Wi (say w) — > (or w — > 1), it can be shown 
that the behaviour of Z is 



ttA(1 + a)A(l + /3)A(-1 - a - (3) 



(5) 



namely, for a (or (3) non negative integer a zero appears, 
while for a, /3 negative integers, we have poles. The 
poles/zeros in a and f3 are connected with behaviour at 
w — > 00. We thus obtain some of the functions A given 
in (0). If we consider the neighbourhood of some of the 



around or 1, 



0, 1, we find that Z behaves as 



-p-p r(l + a + ip) -f-j r(l + [3 + ip) 
Z ~ll T(-a-ip) W' 



T(-P-ip) 



(6) 



where i = 1, 2, • • • , N — 1. 

Using the invariance of Z under w 
also show that Z has the symmetry 



1/w, one can 



Z(a, f3, p) = Z(-a -(3-2-p{N- 1), 0, p) (7) 

As is evident in the canonical partition function (Q) the 
zeros and poles given in terms of p depend on the tem- 
perature. Those zeros and poles are, in general, related 
to phase transitions. It should be noted that at such 
points the free-energy becomes singular. The simplest 
possible phase transition in this system is the collapsing 
of particles into clumps. In order to investigate the col- 
lapsed phase, it is natural to introduce an order param- 
eter describing the average distance of the particles with 
respect to one another. We therefore consider the param- 



eter 



N 

i,j V ij 'l 



which allows investigation of the collapse 



of the system. As it turns out, < r 2 > has several zeros, 
one of them coinciding with a singularity of the parti- 
tion function. We also show that a class of singularites 
of Z may be related to the zeroes of higher moments of 
< r 2 > i.e < r 2q >. Using the invariance of the Hami- 
tonian under translation, we set the coordinate of one of 
the particles to zero and we can write the mean square 
distance as 



< r z >= 



1 < r 2 > 



(8) 



where 



2 



. N-l 

<r 2 > = J ^zn^KI 2 "! 1 "^! 2 ' 3 



N—l 

n H-wj^zf+^i 



z\^\z- Wi \ a p (9) 



Similarly to the case of the partition function, the inte- 
gral displayed in (Q) has a complex structure of poles and 
zeros. Let us classify such a singularity structure, using 
now the hints we have gotten based on the computation 
of the partition function, (||) . 

Supposing w — > (or w —>■ 1, oo) we find 



N-l 



< r 2 > ~ ]J A(l + a + ip)A(l + /3 + i/j) x 

i=l 

A(-l - a - /3 - (AT - 2 + i)p) 



(10) 



At Wat = z — > 0, < r 2 >o behaves as 



T(a+2) 

r(-i-a)- 



For 

i j.— U/ , 

z, w — > 0, it behaves as r^i+a+i/2+p)) • ^ ^ s eas y ^° show 
that analysing the case where several points tend to zero 

we see that < r 2 >o must behave as 



r(l+a+ 7 i T + ip) 
r(l-(a+ 7 i T +ip))' 

where i — 0, 1, • • • , N — 1. The behaviour of < r 2 >o at 
1 similarly obtained. 

Therefore, for the singularity structure of (r 2 ) we arrive 
at the result 

N-2 

(r 2 ) ~ J] A(l + a + ? : P ) 

i=0 

A(l + /3 + ip)A(-l -a-/3-(JV + t- l)p) 

N—l 

Y[ A(l + a + ip- 



1 



t— U 

A(-l - a - /3 - (AT - 1 + i)p - 



1 



N - 1 



:) (11) 



while the normalized 



N—l 



has the form 



1 



<r 2 )~ [] A(l + (i + 2)p- 

i=0 

A(-l- (A^ + 3 + ^)p- 
A(l+(A^-t-l)p) 2 A(-l-2(Ar + l)p) 



i+r 

N-i-l i 



(12) 



for i = it has a zero at 1 + a + 1= positive integer 
= 1,2 - • • therefore the possible value of p is given by 
2p = a = -1,0,1,... . 

The expression for r 2 has several zeros. They can cer- 
tainly not be all physically relevant, since most are just 
consequences of the regularization of the integrals (||). 
After encountering a zero, we loose the physical relevance 
of the integral. Below the critical point, the theory is in 
a clumped phase, while above the particles are far apart. 
Therefore, there should be no further critical point. This 
fixes our critical point to be, at large AT 



kT c = -NGm 2 
4 



(13) 



For large AT this coincides with an old result of Salzberg 
JLH , from wich a possible phase transition point can be 
obtained. 

Note that (r 2 ) behaves as \T — T c \ in terms of the 
temperature. 

The higher moments of < r 2 > have the singularity 
structure 



n 



JV- 
i=0 



-1 A(l+(» + 2)p+ i j T ) 
A(l + (AT + l)p) 



A(-l + (N + 3 + z)p- w ^- I ) 
A(-l - 2(JV + l)p) 



(14) 



There are zeros at the values p = — \q, — \q + 1, • • •• If 
q is a multiple of (i + l), for i = 1, ■ ■ ■ N — 1 there are 
further singularities. Depending on the exact value of q 
the behaviour in terms of the departure from the critical 
temperature changes, which shows that the system has 
multi-fractal nature. The zeros of < r 2q > are coincident 
with a class of singularities of Z. For q = 2(N — 1) the 
behaviour of (r 2q ) is not linear in IT — T r \. 



III. THE GRAND PARTITION FUNCTION 

The grand partition function of the system can be writ- 
ten as 



N=0 ' 



(15) 



where Hjq is given by (||) and z — e M is the fugacity. 

Acoording to [||J^] this many body problem can be 
transformed into a field theoretic one. Using the defini- 
tion of density p(r) — Xa=i H r — it has been shown 
that Zq can be written in terms of Liouville field theory, 



J Dcj)'e T n 



(16) 



where T e // = ^^w~ ano - t 1 ' 2 = z Gm 3 . we rescale <j)' as 
' — ' ~^r't > - Therefore, the action can be written as 



S L = i~ I d 2 x(|V0|^+Me 60 ) 

57T 



(17) 



where b = y 2|1£. an d p = j-^ . In || , it has also been 
shown that the correlations of the density can be written 
in terms of vertex operators of the Liouville theory. 

It is conventional to add term ^-Ry/gip to the Liou- 
ville Lagrangian density, where R is the scalar curvature 
of background metric g M „, and the parameter Q adjusted 
to ensure that all physical quantities be independent of 
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a particular choice of the background. However it is pos- 
sible to choose a specific background which is flat every- 
where except for few selected points |H| . 

The Liouville field <fi(z, z) is a logarithmic operator fl(f | 
and varies under holomorphic coordinate transformation 



z — > w(z) as 4>{w,w) 



(z,z) 



log 



dw 12 

dz I 



where Q 



^ + i ■ Q parametrizes the central charge of the theory 
by the well knwn relation c = 1 + 6Q 2 . 

The introduction of Q in the Liouville action is done 
for taking into account the fact that the theory corre- 
sponds to a Coulomb gas where there is no zero charge 
sector, except for the vacuum. It corresponds to differ- 
ent boundary conditions, and takes account of the zero 
modes of the theory |8j . In the description of ref . || there 
is no zero mode in the inverse propagator, fixed as being 
the classical gravity potential, which naturally matches 
the fact that the background curvature has to vanish in 
Minkovski space. 

We suppose that the theory can be described as a com- 
pactified Euclidian space, in which case the Q-tevm. can 
be seen as a requirement of renormalization. The ques- 
tion to be answered is whether the large compactification 
radius limit is smooth or not. Since our argumentation 
are based on local Green functions, as in (|2^) below, we 
do not expect any major difference to occur. On the 
other hand, the argumentation based on the conformal 
dimensions (as in (24)) below and following equations) 
may depend on the above limit. However, it is rewarding 
to see that the results are in accordance with general ex- 
pectations, leading to the conjecture that the large com- 
pactification radius limit is smooth. 

It is nevertheless necessary to stress that the introduc- 
tion of the Q term is non trivial, and might lead to a 
change of the problem. The aim here is to describe the 
results obtained comparing them with known results. We 
know in fact from || that Liouville theory describes grav- 
ity, and any effort in the direction of understanding the 
model is worth undertaking. 

It is well known that the exponential Liouville opera- 
tors V a (x) = e 2av ^ are the spinless primary conformal 
fields of dimension A Q = a(Q — a). The two, three and 
four-point correlation functions of Liouville field theory 
for given a have been calculated in ||[l5| . 

In addition, we find that the dependence on the scale 
p, of any correlation function in Liouville theory [|l7 11 
is 



(tt/z) 



(Q-EjLia 4 )/6 



(18) 



We note that the Liouville theory is self-dual under 
b — ► 1/6. Indeed, we consider the partition function of the 
Liouville theory. According to || the partition function 
has the form 



Z = - 



7T/HX& 2 ) . V 6" r(-i/6 2 ) 



V2ir 2 (b+l/b) y T{l-b 2 )> r(l/6 2 -l) 



(19) 



Under transformation of b and p as given by 



i 



ttA(1/6 2 



-(^7(n) 



2NNl/b 2 



(20) 



it is easy to show that the partition function transforms 
as 



(21) 



This duality transformation was first observed by 
Zamolodchikov jl5| . We further observe that there exists 
a sequence of critical values for 6, so that the partition 
function becomes singular, that is 



N 



(22) 



where b c = 1. As in the case of the canonical partition 
function, the grand-canonical Zq has zeros. 

In order to understand the nature of phase transition 
in the grand-canonical ensemble we define the variance 
of the number of particles as a parameter, where in the 
clumping point (AA) 2 = (< A 2 > - < A > 2 ). This 
variance has a peculiar behaviour. It is easy to show that 
the (A A) 2 can be witten in terms of two-point density 
correlation functions as: 



(AA) 2 



d 2 xid 2 x 2 < p(xi)p(x 2 ) > 



(23) 



Now using the results of || for two point correlation func- 
tions of Liouville vertex operators, one can show that at 
the transition point (p(xi)p(x2)) —* 0, which shows that 
particles collapse to local clusters. 

Now it is possible to investigate the multi-fractal na- 
ture of 2D-Gravitating gas considering moments of the 
density in a given scale R. Using the conformal dimen- 
sions of Liouville operators we can show that 



d 2 x < p q >~ R- T(q) 



(24) 



where r(q) = 2(1 — b 2 q)(q — 1) which is valid only for 



< 



1-2 

Equation ( |24| ) allows us to determine the distribution 
functions of density p(p) , such that 



and 



(P n ) = / P n p(p)d P 
Jo 



p(p) = f(p)e~^^ log {pR } 



(25) 



(26) 



where f(p) is a smooth function of log(p). Moreover 
a = 166 6 , d = 26(1 + 46 2 ) (see also Q). This is the 
famous log-normal distribution, considered as a charac- 
teristic feature of a disordered system. It turns out that 
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for a dilute gas p ~ ^, 26(1 1 +4b y ) , and we find the Gaussian 
distribution functions, whose variance is controlled by R. 
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